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Pairing Effects on the Fragment Mass Distribution of Th, U, Pu, and Cm Isotopes* 
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' Department of Physics, Liaoning Normal University, Dalian 116029, China 


This article presents a comprehensive study on the fission process in Th, U, Pu, and Cm isotopes using a 
Yukawa-folded mean-field plus standard pairing model. The focus is on analyzing the effects of the pairing 
interaction on the fragment mass distribution and its dependence on nuclear elongation. The study demon- 
strates that the pairing interaction plays an important role in the fragment mass distribution for ??°Th, 73*U, 
240Pu, and *46Cm. Numerical analysis reveals that increasing the pairing interaction strength leads to a de- 
crease in asymmetric fragment mass distribution and an increase in symmetric distribution. Furthermore, the 
investigation examines the odd-even mass differences at symmetric and asymmetric fission points, highlighting 
their sensitivity to changes in the pairing interaction strength. The systematic analysis of Th, U, Pu, and Cm 
isotopes’ fragment mass distributions demonstrates the effectiveness of the model in reproducing experimental 
data. Additionally, the study explores the effects of the zero-point energy parameter and neck-breaking prob- 
ability parameter on the fragment mass distribution for 24°Pu. In conclusion, this research provides valuable 
insights into the fission process by emphasizing the importance of the pairing interaction and its relationship 


with nuclear elongation. 
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I. INTRODUCTION 


Nuclear fission is a fundamental process that plays a central 
role in modern nuclear technology. The theoretical calcula- 
tion of the fission process is a complex and challenging prob- 
lem, necessitating the utilization of advanced nuclear models 
and computational techniques [1—6]. Throughout the years, 
numerous theoretical models have been developed to predict 
fission yields, ranging from simple empirical models to so- 
phisticated microscopic models based on nuclear structure 
and reaction theory [7-9]. These models have been validated 
against experimental data and proven to be valuable tools for 
predicting the behavior of nuclear systems. 

Pairing interactions have a significant impact on the prop- 
erties of the fissioning nucleus and the resulting fission prod- 
ucts. For instance, the strength of the pairing interaction 
strongly influences the shape of the barriers separating the 
ground state from scission [10-21] fission fragment distri- 
butions [22-25], and spontaneous fission lifetimes [26]. In 
the dynamical description of nuclear fission, pairing interac- 
tion should be considered on the same footing as those as- 
sociated with shape degrees of freedom [15]. Understanding 
the role of pairing interactions in nuclear fission is an active 
area of research, and various theoretical models have been de- 
veloped to describe their behavior in different fission scenar- 
ios. Macroscopic-microscopic studies have demonstrated that 
pairing fluctuations can significantly reduce collective action 
and affect the predicted spontaneous fission lifetimes [27]. 
In the Hartree-Fock-Bogoliubov (HFB) model, pairing can 
be self-consistently included by extending the trial space to 
quasi-particle Slater determinants [22, 28]. Theoretical works 
based on the HFB method have revealed that the effect of 
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pairing interactions hinders collective rotation, reduces level 
crossings, and thus shortens the half-lives of spontaneous fis- 
sion [29]. The role of dynamical pairing in induced fission 
dynamics has been investigated using the time-dependent 
generator coordinate method in the Gaussian overlap approx- 
imation, based on the microscopic framework of nuclear en- 
ergy density functionals [30]. It has been shown that the in- 
clusion of dynamical pairing has a pronounced effect on col- 
lective inertia, the collective flux through the scission hyper- 
surface, and the resulting fission yields. The latest research 
on the fission dynamics mechanism of 74°Pu, based on the 
time-dependent Hartree-Fock method, demonstrates that as 
dynamical pairing diminishes at high excitations, the random 
transition between single-particle levels around the Fermi sur- 
face, to mimic thermal fluctuations, becomes indispensable in 
driving fission [31]. 


In recent studies, an iterative algorithm [32, 33] has been 
employed to investigate the fission barriers and static fission 
paths of Th, U, and Pu isotopes using the deformed mean- 
field plus standard pairing model with an exact pairing solu- 
tion [34]. This innovative approach provides a precise rep- 
resentation of pairing interactions in nuclear fission, avoid- 
ing artifacts introduced by BCS calculations, such as non- 
conservation of particle numbers and pairing collapse phe- 
nomena [11]. The comprehensive investigation of the inner 
and outer fission barriers in even-even nuclei Th, U, and Pu 
isotopes clearly demonstrates the standard pairing model’s 
ability to closely replicate the experimental inner and outer 
barrier heights in comparison to the BCS scheme [34]. More- 
over, employing the deformed mean-field plus standard pair- 
ing model, researchers have explored the influence of pair- 
ing interactions on scission configurations, total kinetic en- 
ergy, and mass distributions of U isotopes [35]. The model 
successfully reproduces the total kinetic energy and fragment 
mass distributions of ?3?—238U isotopes, exhibiting excellent 
agreement with experimental data. The results highlight the 
scission region’s sensitivity to variations in pairing interac- 
tion strength, particularly for asymmetric and symmetric scis- 


sion points. Notably, changes in the peak-to-valley ratio of 
the mass distribution, resulting from variations in the pair- 
ing interaction strength, underscore the significant impact of 
pairing interactions on the fission process of 2°°U within this 
model. 

It is crucial to develop a reliable and valid model to char- 
acterize the fragment mass distribution, and actinide nuclei 
serve as important regions to test the reliability of these mod- 
els. Therefore, extending our previous work to describe ac- 
tinide nuclei and studying the influence of interactions on the 
fragment mass distribution is both necessary and meaningful. 

This investigation presents a systematic analysis of fission 
fragment mass distributions in Th, U, Pu, and Cm isotopes 
using the deformed mean-field plus standard pairing model. 
The potential energy is calculated within the macroscopic- 
microscopic framework, incorporating the Fourier shape 
parametrization combined with the LSD model + Yukawa- 
folded potential. The mass distribution of fission fragments 
is described using the three-dimensional collective model of 
the Born-Oppenheimer approximation (BOA). Building upon 
our previous work in Ref. [35], this study offers a comprehen- 
sive analysis of the pairing’s impact on the mass distribution 
of fission fragments across the Th, U, Pu, and Cm isotopes 
chain. 


Il. THEORETICAL FRAMEWORK AND NUMERICAL 


DETAILS 
A. Deformed mean-field plus standard pairing model 


The Hamiltonian of the deformed mean-field plus standard 
pairing model for either the proton or the neutron sector is 
given by 


() 


100 


H= yah -GX S$ Sp. 
4=1 it! 


Here, the sums run over all given i-double degeneracy levels 
of total number n, G' > 0 represents the overall pairing inter- 
action strength. The single-particle energies ei are obtained 
from mean-field methods, such as Yukawa-folded single- 
particle potential, Woods-Saxon potential (WS), Hartree- 
108 Fock (HF). The fermion number operator for the i-th dou- 
107 ble degeneracy level is defined as n; = al aip + al a; j 
and the pair creation (annihilation) operator is represented by 
SF = alal, [S; = (SF)! = aj ait]. The up and down 
arrows in these expressions refer to time-reversed states. 

Using the Richardson-Gaudin method [36—41], the exact 
k-pair eigenstates of (1) with v; = 0 for even systems or 
113 Vy = 1 for odd systems, where i’ is the label of the double 
114 degeneracy level occupied by an unpaired single particle, can 
115 be expressed as: 
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117 where |v; ) is the pairing vacuum state with the seniority v; 
18 that satisfies S; |v) = 0 and Ts |i) = Oi Vi |Vir) for all i. 


19 Here, € is an additional quantum number for distinguishing 


120 different eigenvectors with the same quantum number k and 
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1 
121 st (x) = 5 EOT oei (3) 
jai Ue — 26; 
122 in which the spectral parameters ai’) (u=1,2,...,k) satisfy 
123 the following set of Bethe ansatz equations (BAEs): 
oF ii 1 
m 1+0} o -26 YO ga Eia 
é Tu — 2€; u'=1 (£u) Tu = Tu 


where the first sum runs over all ¿ levels and Q; = 1 — ĝi vi. 
For each solution, the corresponding eigenenergy is given by 


k 
B® = y a) + vey. (5) 
p=1 


The general method to find solutions of Eq.(4) is based 
on the polynomial approach described in Refs.[42—45]. This 
approach involves solving the second-order Fuchsian equa- 
tion [46], given by: 


A(x)P"(x) + B(x)P'(x) — V(x) P(x) (6) 


(xi) — 2e;) is an n-degree polynomial, 


where A(x) = JĮ: 


i=l 


Qi 1 
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B(z)/Als) =— Dm (7) 
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The polynomials V(x), also known as Van Vleck polyno- 
mials [46], are of degree n — 1 and are determined based on 
137 Eq. (6). They are defined as follows: 


(8) 


138 


The polynomials P(x) with zeros corresponding to the so- 
lutions of Eq. (4) is defined as 


D 
a] 


(9) 


Here, k represents the number of pairs, and b; and a; are 
the expansion coefficients that need to be determined instead 
of the Richardson variables x;. Additionally, when we set 
ak = 1 in P(x), the coefficient ag—1ı is equal to the negative 
sum of the P(x) zeros, i.e., @aķ-1 = — > ao) = ae! 

For doubly degenerate systems with Q; = 1, if the value of 
x approaches twice the single-particle energy of a given level 
6, i.e., x = 2€5, we can rewrite Eq.(6) as follows[42, 45]: 
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The iterative algorithm for obtaining the exact solution of 
the standard pairing problem using the Richardson-Gaudin 
method is established, employing the polynomial approach 
described in Eq.(10)[32]. This algorithm showcases remark- 
able efficiency and robustness, capable of handling both 
spherical and deformed systems on a large scale. A crucial 
element contributing to its success lies in the determination 
of initial guesses for the large-set nonlinear equations, en- 
suring control and adherence to fundamental physical prin- 
ciples. Moreover, the algorithm effectively tackles the chal- 
lenges of nonsolutions and numerical instabilities frequently 
encountered in existing approaches by reducing the high- 
dimensional problem to a one-dimensional Monte Carlo sam- 
pling procedure. Leveraging this innovative iterative algo- 
rithm, we employed the model to explore actinide nuclei iso- 
topes, resulting in exceptional agreement with experimental 
data [32-35]. 


B. The Fourier shape parametrization 


Recent studies have highlighted the remarkable efficiency 
of the Fourier parametrization when describing the essential 
features of deformed nuclear shapes, extending up to the scis- 
sion configuration [7, 47]. Building upon these findings, the 
present work employs the innovative Fourier parametrization 
of nuclear shapes in conjunction with the LSD + Yukawa- 
folded macroscopic-microscopic potential-energy prescrip- 
tion, yielding highly efficient results [35, 48, 49]. Specifi- 
cally, the macroscopic-microscopic framework introduced in 
Ref. [35] serves as the primary foundation for this study. In 
this framework, the single-particle energies £; in the model 
Hamiltonian (1) are derived from the Yukawa-folded poten- 
tial. The expansion of the nuclear surface, expressed as a 
Fourier series in terms of dimensionless coordinates, is given 
by: 


p(z) 
Re 


(2n — 1)r z — a 


ji 


where p°s(z) represents the distance from a surface point to 
the symmetry z-axis, and RO = 1.2A!/3 fm corresponds to 
the radius of the corresponding spherical shape with the same 
volume. The shape extends along the symmetry axis by 220, 
with the left and right ends located at zmin = Zsh — Zo and 
Zmax = Zsh +20, respectively. Here, zo is half of the extension 
of the shape along the symmetry axis, derived from volume 
conservation, and Zs is determined to ensure the center of 
mass of the nuclear shape is positioned at the origin of the 
coordinate system. Following the convergence properties dis- 
cussed in Ref. [7], we retain the first five orders az,..., ag 
as a Starting point, and transform the parameters an into the 


= > [azn cos ( 
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deformation parameters qn as follows: 
q= a9” /az — az/a8”, 
q3= 43, 
q4= a4 + V (02/9)? T (a[)2, 
g5= a5 — (qa = 2)a3/10, 
q6= 46 — V (q/100)? + (ag”)2, 


(12) 
where a? representing the values of the Fourier coefficients 
for the spherical shape. While the higher-order coordinates q5 
and qe are typically negligible within the current approach’s 
accuracy, the set of q; possesses explicit physical significance 
in characterizing the nuclear fission process. Specifically, q2 
denotes the elongation of the nucleus, q4 represents the neck 
parameter, and q3 reflects the left-right asymmetry parame- 
ter. For this study, the dynamic process of nuclear fission 
will be modeled in the three-dimensional deformation space 
(q2, 93; q4) using the Fourier shape parametrization. It is im- 
portant to note that the present work does not consider nonax- 
ially symmetric shapes since they primarily play a significant 
role in the vicinity of the ground state and the first saddle 
point. 


C. The mass distributions 


In previous studies, the use of Wigner functions to ap- 
proximate the probability distribution associated with neck 
and mass asymmetry degrees of freedom has shown good 
agreement between model predictions and experimental re- 
sults [7, 49-52]. Building on these ideas, this study proposes 
a fission dynamics scenario where the motion towards fission 
primarily occurs along the q2 direction, accompanied by fast 
vibrations in the perpendicular q3 and q4 collective variables. 
The total eigenfunction Vn g (q2, q3, q4) of the fissioning nu- 
cleus is then approximated as the product of two functions: 


Une (92,93, 94) = Hne (92) On (93; q4; 92). (13) 


Line (q2), which depends mainly on a single variable q2 
and describes the motion towards fission, and ¢,,(q3, q4; q2), 
which simulates n-phonon fast collective vibrations on the 
perpendicular 2D plane q3, q4 for a given elongation q2. For 
low-energy fission, only the lowest energy eigenstate @n=o is 
considered. 

The density of probability W(q3,q1;q2) of finding the 
system for a given elongation q2 ,within the area of (q3 + 
dq3, q4 + dq4), is given as 


W (q3, q4; 92) = [Y (a2, 93, 94)|? = |0(a3, qa; G2)? 14) 


To take into account the fission process, a Wigner function 
is employed, which is given by: 


V(q3, q4; q2) — Vmin(q2) 
Eo 


W (q3, q4; 2) & exp{ kas) 


244 Where Vmin(q2) is the minimum potential for a given elon- 
245 gation q2, and Eo is the zero-point energy, treated as an ad- 
24 justable parameter. 

To obtain the fragment mass yield for a given elongation 
248 q2, the probabilities coming from different neck shapes, sim- 
249 ulated by the q4 parameter, are integrated as 


247 


250 


w(q3; 42) = [wo qa; q2)dqa. (16) 
251 Following the concept introduced in Ref. [51], the neck rup- 
ture probability P is assumed to be equal to 

ko 
253 P(q3, q4, 92) = 7 Freck (Rneck), (17) 
where k represents the momentum in the direction towards 
fission, and the constant parameter ko serves as a scaling pa- 
256 rameter. Ryeck is the deformation-dependent neck radius, 
and Preck iS a geometrical factor indicating the probabil- 
ity of neck rupture, which is proportional to the neck thick- 
ness. The expression for the geometrical probability factor 
Preck(Rneck) can be chosen arbitrarily to some extent, such 
as using Fermi, Lorentz, or Gaussian functions [52]. In this 
study, we adopt the Gaussian form as follows: 


Preck(Rneck) = exp[—In2(Rneck/d)”], 


where d represents the half-width of the probability distribu- 
tion, and it is treated as another adjustable parameter in this 
analysis. The momentum k in Eq. (17) simulates the dynam- 
ics of the fission process, which, as usual, depends on both 
the local collective kinetic energy E — V (q2) and the inertia 
towards the leading variable q2. 

re 


2M (q2) 


(18) 


= Ekin = E —Q—V (Q2), (19) 


where M (q2) represents the averaged inertia over the de- 
grees of freedom q3 and q4 at a given elongation q2, while 
V (q2) denotes the averaged potential. Here, it is assumed 
that the portion of the total energy converted into heat, de- 
noted as Q, is negligibly small. A convenient approximation 
for the inertia M (q2) is to employ the irrotational flow mass 
277 parameter Birr [53], which is initially expressed as a func- 
tion of the single collective parameter R12, representing the 
distance between fragments, and the reduced mass u of both 
fragments. 


— ƏR \" 
M (q2) = mel + 11.5(Birr/p = 1)] a . (20) 
q2 
To incorporate the neck rupture probability P(q3,q4;q2) 
in Eq.(17), the integral over the probability distribution in 
Eq.(15) with respect to q4 needs to be reformulated. This is 
achieved by expressing it as follows: 


w(q3; q2) = f Wlas a 42) Plaz a, as)das (21) 


The aforementioned approximation implies a crucial obser- 
vation: for a fixed q3 value, fission may occur within a spe- 
cific range of q2 deformations, each associated with different 


290 probabilities. To obtain the accurate fission probability distri- 
bution w’ (q3; q2) at a particular q2 value, it is necessary to ex- 
clude fission events that occurred in previous configurations 
with q4 < q2. This can be achieved by applying the following 
expression: 


291 


= Sys qq (935 1) AA 
f was; dg 


w’ (93; q2) = w(q3; q2) (22) 


The normalized mass yield is then obtained as the sum of 
297 partial yields at different given qo: 


J w (933 2) dae 
Y (q3) = , 
(a) J w' (q3; q2)dq3dq2 


(23) 


Since the scaling parameter ko introduced in Eq. (17) does not 
longer appear in the definition of the mass yield, the only free 
parameters, the zero-point energy parameter Eo in Eq. (14) 
and the half-width parameter d appear in the probability of 
neck rupture (18). The free parameters in the model are the 
zero-point energy parameter Eo and the half-width parame- 
ter d, with d = 0.16Rp and Ey = 2.2 MeV used in this 
work, based on their successful reproduction of experimental 
fragment mass yields in low-energy fission of Pt to Ra iso- 
topes [49]. 


Il. THE POTENTIAL ENERGY 


In this study, the potential energy of the system is com- 
puted using the macroscopic-microscopic approach. The total 
energy of a nucleus with a specific deformation, represented 
as Erota (N, Z, dn), is determined through the following pro- 
cedure: 


Erota (N, Z, dn) = Exp(N, Z) ale Ep(N, Z, Qn): (24) 


In the calculation, the total energy Erotaı( N, Z, dn) is com- 
317 posed of two main contributions. The first term, denoted 
as ELp(N, Z), corresponds to the macroscopic energy cal- 
culated using the standard liquid drop model. This term 
takes into account the proton number Z and neutron num- 
ber N [54]. The second term, Ep(N, Z, qn), is related to 
the shape parameters gz, q3, qa and represents the potential- 
energy surface. In the current calculation, the focus is solely 
on this energy term, neglecting other contributions to the total 
energy. 


Ep(N, Z, qn) = Eaet (N, Z, Qn) T Esen (N, Z, Qn): 


ae Epair(N, Z, Qn) (25) 


The deformation correction energy Euer(N, Z, q2, 93, 94) 
is obtained from tables in Ref. [55]. The micro- 
scopic terms consist of the shell correction energy 
aa (N, Z, {ei}, q2, 93, q4) proposed by Strutin- 
sky [56, 57] and the pairing interaction energy 
333 a (N, Z, {€:},92,93,q4) calculated using Eq. (1), 
334 Where v (m) represents the label of the neutron (proton) 
335 Sector. The microscopic calculations involve considering 18 
sss deformed harmonic oscillator shells in the Yukawa-folded 
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Fig. 1. (Color online) Contour map of the potential-energy surface 
of the nucleus 74°Pu (in MeV), minimized q4 with the pairing inter- 
action strength G” = 0.08 and G7 = 0.10 (in MeV). The black 
trajectory shows the static fission path. 


single-particle potential to determine the single-particle 
energy levels. Additionally, for the pairing correction energy, 
66 single-particle levels around the neutron Fermi level 
and 51 single-particle levels around the proton Fermi level 
are taken into account. To determine the overall potential- 
energy surface, a multidimensional minimization process 
is performed, simultaneously considering all axial degrees 
of freedom. This includes minimizing the elongation of 
the nucleus q2, the asymmetry of the left and right mass 
fragments q3, and the size of the neck q4. By considering 
all these degrees of freedom together, a comprehensive 
understanding of the nuclear shape and energy landscape can 
be obtained. 

Figure 1 illustrates the behavior of the potential-energy 
surface (PES) during the fission process for 74°Pu. At the ini- 
tial stage of fission q2 < 0.5, the the PES exhibits a very soft 
octupole deformation, and the minimum of the PES (ground- 
state) is located at q) = 0. The fission barrier heights ob- 
tained from the present model are in good agreement with the 
corresponding experimental results from Ref. [58]. Specif- 
ically, the inner barrier height is 4.88 MeV, the outer bar- 
rier height is 5.24 MeV, while the experimental results are 
5.80 MeV and 5.30 MeV, respectively. Furthermore, in the 
asymmetric fission path, Figure 1 shows a plateau at high de- 
formation followed by a cliff (the asymmetric scission point: 
q2 = 2.45, q3 = 0.10, q4 = —0.09). 

The strength of the pairing interaction G is typically de- 
termined using empirical formulas or by fitting experimental 
data, such as the odd-even mass difference [59-62]. Pre- 
vious studies have demonstrated that pairing plays a crucial 
role in the region of the inner and outer barriers, and that 
the first and second saddle points are highly sensitive to the 
strength of the pairing interaction [34, 64, 65]. Therefore, in 
the present model, experimental observables such as the odd- 
even mass difference (reflecting ground-state properties) and 
barrier heights (reflecting excited-state properties) are used to 
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Table 1. Pairing interaction strength G” (G”) (in MeV) for Th, U, 
Pu, and Cm isotopes. 


Th U Pu Cm 
G” 0.096 0.080 0.080 0.096 
G™ 0.120 0.100 0.100 0.120 
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Fig. 2. (Color online) Odd-even mass differences (in MeV) for Th, 
U, Pu, and Cm isotopes. Experimental values are denoted as "Expt." 
and the theoretical values calculated in the present model are de- 
noted as "Theor.". Experimental data are taken from Ref. [66] (in 
MeV). 


determine the experimental values of the pairing interaction 
strength in the fission process. 

In this paper, realistic values of the pairing interaction 
strengths for the isotopic chains of Th, U, Pu, and Cm were 
obtained by fitting the experimental values of the odd-even 
mass difference and the heights of the inner and outer bar- 
riers. The odd-even mass difference is calculated using the 
three-point formula: 


P(A) = Erota (N ac 1, Z) T Erota (N = 1, Z) 

— 2Etota(N, Z). (26) 
The odd-even mass difference is attributed to the presence 
of nucleonic pairing interactions and is highly sensitive to 
changes in the pairing interaction strength G [66]. The corre- 
sponding values of G” (G7) are shown in Table 1. 

Figure 2 clearly demonstrates that the odd-even mass dif- 
ferences obtained using the present approach closely match 
the experimental data for Th, U, Pu, and Cm isotopes. Ad- 
ditionally, as depicted in Fig. 3, the inner fission barriers (a) 
for Th, U, Pu, and Cm isotopes, as well as the outer fission 
barriers (b) for the same isotopes, calculated in the current 
model exhibit remarkable agreement with the corresponding 
experimental values. It is necessary to indicate that the theo- 
retical inner barrier heights of light Th isotopes in Fig. 3-(a) 
are systematically lower than the experimental data, which 
is also reported in other calculations for light actinides in 
Refs. [12, 13, 67-69]. Based on analysis of different effects of 
neutron and proton pairing interaction on the inner and outer 
barrier heights in Ref. [32], the result above may be related to 
the strong neutron pairing interaction strength. In this article, 
the pairing interaction strength value in Table 1 is set to GG 
(GG) for Th, U, Pu, and Cm isotopes. 
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Fig. 3. (Color online) Inner fission barriers (a) and outer fission bar- 
riers (b) for Th, U, Pu, and Cm isotopes, respectively. Theoreti- 
cal values obtained using the present model are labeled as "Theor. ", 
while experimental values are indicated as "Expt.". The experimen- 
tal data is sourced from Refs.[12] and is measured in MeV. It is im- 
portant to note that the typical uncertainty in the experimental val- 
ues, estimated based on variations among different compilations, is 
approximately +0.5 MeV[12]. 


IV. EFFECT OF PAIRING INTERACTION ON THE 
FRAGMENT MASS DISTRIBUTIONS OF ?°°Th, ?4 U, 
240 Pu , AND °? Cm 


The investigation of dynamics around fission structures 
plays a crucial role in comprehending various aspects of the 
final fission state, such as kinetic energy and mass distribu- 
tions [7, 70, 71]. In this study, the fission fragment mass dis- 
tribution of ?40Pu was calculated based on its potential energy 
surface and compared with experimental data [72]. 

As shown in Fig. 4, the agreement between the calculated 
results and experimental data is reasonable [72]. Moreover, 
the obtained fission fragment mass distribution aligns with 
the understanding that static fission in 74°Pu is predominantly 
asymmetric, as indicated by the fission potential energy sur- 
face. In the calculation process, a Gaussian folded function 
with a full width at half maximum (FWHM) of 4.9u [73] is 
employed to determine the mass yields. Additionally, two pa- 
rameters, Eo = 2.2 MeV and d = 1.6fm [49], are utilized. 

In the case of fission nuclei, each elongation deformation 
variable q2 corresponds to a distribution of fragment mass 
numbers Ap of the nuclear fragments produced during fis- 
sion. Figure 5 illustrates the distribution of fragment mass 
numbers for 74°Pu. It is observed that fission predominantly 
occurs in the region of asymmetric fission, with the cor- 
responding mass numbers of the heavy fragments centered 
around A ~ 141. The scission point, representing the point 
of fragment separation, is located at q2 = 2.3. Only a small 
portion of fragments undergo symmetric fission. 

To investigate the influence of the pairing interaction on the 
fission fragment mass distribution under the current model, 
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Fig. 4. (Color online) Mass yield of 7“°Pu and compared with ex- 
perimental data [72]. 
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Fig. 5. (Color online) Mass yield of 7*°Pu as a function of the mass 
number Ay, and the elongation deformation q2. 


we calculated the yield of the fission fragment as a function of 
mass numbers (Ap) for 230TH, 2340, 240Pu, and 246Cm with 
different pairing interaction strengths. The results depicted in 
Fig. 6 indicate that for these nuclei, as the pairing interaction 
strength G increases from 80%G‘ to 120%Go, the two asym- 
metric peaks of the theoretical yield are significantly reduced, 
while the symmetric valley becomes more prominent. Similar 
observations were reported in a three-dimensional Langevin 
model based on the BCS approximation [74]. These find- 
ings suggest that the fragment mass distribution is sensitive to 
variations in the pairing interaction strength and highlight the 
significant role of pairing interaction in determining the frag- 
ment mass distribution for 2°°Th, 224U, 24°Pu, and 24°Cm. 
Furthermore, when the pairing interaction strength G is set 
to 120%Gb, the theoretical calculations closely match the ex- 
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perimental data for the fragment mass distribution of 7°°Th, 
234U, and 249Pu. However, for 24°Cm, the calculated results 
align better with the experimental values at a pairing interac- 
tion strength of 80%Go. 
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Fig. 6. (Color online) Mass yields for 230Th, ?34U, ?40Pu, and 
246 Cm as a function of mass numbers (Af) with varying pairing in- 
teraction strengths. Experimental data for ?°°Th are extracted from 
the charge-yields as reported in Ref [65]. The mass yields for ?°4U 
are obtained from Ref [75]. For 24°Pu , the calculated mass yields 
are compared with experimental data [72]. The experimental data 
for 24°Cm are taken from Ref. [76]. 


Figure 7 illustrates the calculated odd-even mass differ- 
ences at the asymmetric and symmetric fission points for 
230Th, 234U, 240Pu, and 246Cm, considering the variation of 
the pairing strength G ranging from 80%G to 120%Go. In 
this analysis, it is assumed that the ground-state odd-even 
mass differences represent the odd-even binding-energy dif- 
ferences in the scission configuration, despite some shape 
differences. By comparing the experimental odd-even mass 
differences of asymmetric and symmetric fission fragments 
in nuclei such as 28°Th, ?34U, 24°Pu, and 246Cm, the calcu- 
lated results show better agreement with experimental values 
at the asymmetric fission point when the pairing strength is 
set to 120%Go for ?°°Th, 734U, and 74°Pu. Conversely, at 
the symmetric fission point, a stronger pairing interaction is 
required, and the calculated results agree better with exper- 
imental values when the pairing strength is set to 140%Go 
for 72°Th, ?34U, and 249Pu. The calculated results for the 
odd-even mass differences at the symmetric and asymmetric 
fission points for ?46Cm demonstrate that pairing strengths of 
80%Go and 120%Go match well with experimental values. 
This finding aligns with the earlier conclusion that increasing 
the pairing interaction strength leads to a better distribution 
of fission fragment masses for for 7°°Th, ?34U, and 74°Pu. 

The calculations presented above suggest that to achieve 
a better description of fission products, different elongation 
deformations of nuclei require different strengths of the pair- 
ing interaction. By fitting the ground-state binding energy, 
inner and outer barrier heights, and mass distribution calcu- 
lations for Pu isotopes, the optimal values for the strength of 
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Fig. 7. (Color online) Odd-even mass differences (in MeV) of 


230Th, ciel OF 240 py, and 74°Cm at the asymmetric scission point 
and the symmetric scission point with varying pairing interaction 
strengths G’“*) from 80%Go to 120%Go (in MeV). The theoret- 
ical values calculated in the present model based on Eq. (26) in 
Ref. [60] are represented as "Sym.Theor." for symmetric point and 
"Asym.Theor." for asymmetric point. The experimental values of 
the odd-even mass difference for asymmetric and symmetric fission 
fragments of 230TH, 234U, 240Pu, and 24°Cm denoted as "Expt." are 
obtained from Ref. [66] (in MeV). 


ass the pairing interactions have been determined. As depicted in 
486 Fig. 8, the strength of the pairing interactions exhibits a non- 
487 linear variation with increasing elongation deformation of the 
488 nucleus. In comparison to the barrier height, a stronger in- 
4s9 teraction is required to accurately describe the fragment mass 
490 distribution. 
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Fig. 8. (Color online) Pairing interaction strength G” (G7) (in 
MeV) obtained by fitting the ground-state binding energy, inner and 
outer barrier heights, and fragment mass distribution calculations for 
236—242 Pu isotopes. Points (A)-(E) represent the corresponding q2 
values for the ground-state binding energy, inner and outer barrier 
heights, and the asymmetric and symmetric scission points. 
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V. FRAGMENT MASS DISTRIBUTION OF Th, U, Pu, AND 
Cm ISOTOPES 


Based on the above results, this work calculates the frag- 
ment mass distribution of Th, U, Pu, and Cm isotopes chain 
based on the corresponding potential energy surface, with the 
pairing interaction strength set to 120%Go. The theoreti- 
cal calculations, as shown in Fig. 9, exhibit good agreement 
with experimental data for all isotopes. The peak height, peak 
width, and peak position of the fragment mass distribution 
closely match the experimental data. However, there are some 
discrepancies for specific isotopes, which can be attributed to 
the limitations of available experimental data. 
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Fig. 9. (Color online) Mass yields for Th, U, Pu, and Cm isotopes as 
a function of mass numbers (A+). Theoretical values calculated us- 
ing the present model are represented as "Theor." while experimental 
data is denoted as "Expt." [75]. The experimental data for ??°Th and 
230Th are obtained by converting the charge distribution with an ex- 
citation energy of 11 MeV [65]. For the isotopes ?°U, ?°°U, and 
242 Pu, the experimental data used is from thermal neutron-induced 
fission [77], while for 236Pu, 238Pu, and 74°Pu, the data is from 
spontaneous fission experiments [76]. The evaluated post-neutron 
data for 7°°U and 744-24°Cm are taken from ENDF/B-VIII.O [76]. 


For 725Th and 73°Th, the experimental data for the frag- 
ment mass distribution were obtained by converting the 
charge distribution of the fragments at an excitation energy 
of 11 MeV in the fission system [65]. This may explain why 
the experimental value of the asymmetric mass yield of ??°Th 
is lower than the theoretical value, while the symmetric fis- 
sion yield is relatively high. Experimental data from thermal 
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neutron-induced fission were used for ?34U, ?36U [75]. The 
theoretical results show a higher symmetric valley for ?34U 
compared to the experimental data. For 7°°U, due to a lack of 
available experimental data, evaluated post-neutron data from 
ENDF/B-VIUII.0 were utilized [76]. 

Experimental data from spontaneous fission were used for 
236Pu, 238Pu, and 74°Pu, and the calculated results closely 
match the experimental data in terms of peak width. For 
242Pu, experimental data from thermal neutron-induced fis- 
sion were employed, and the calculated results exhibit similar 
peak width but deviate from the experimental data by 2 — 3 
mass units in peak position. For 244~?48Cm isotopes, evalu- 
ated post-neutron data from ENDF/B-VIII.0 were used. The 
calculated results in Fig. 9 demonstrate good agreement with 
experimental data, indicating the effectiveness of the present 
model in reproducing the fission fragment mass distribution. 

Overall, the model employed in this work successfully re- 
produces the experimental data of the fission fragment mass 
distribution for Th, U, Pu, and Cm isotopes, providing a valu- 
able tool for understanding and analyzing fission processes. 


VI. EFFECTS OF MODEL PARAMETERS ON 
FRAGMENT MASS DISTRIBUTION OF 2*°Pu 


In the subsequent research, the effects of the zero-point en- 
ergy Eo in Eq. (15) and the neck-breaking probability param- 
eter d in Eq. (18) of the three-dimensional collective model 
on the fragment mass distribution of 74°Pu were investigated. 
The results, as depicted in Fig. 10-(a), indicate that the 
neck-breaking probability parameter d primarily influences 
the peak position of the asymmetric peak. When the neck- 
breaking probability parameter d increases, the peak posi- 
tion of the asymmetric peak shifts towards larger fragment 
masses. 

On the other hand, the zero-point energy Ey mainly af- 
fects the peak value of the fission fragments. As depicted 
in Fig. 10-(b), when the zero-point energy Eo increases, the 
asymmetric peak value of the fission fragment mass distri- 
bution decreases. These observations are consistent with the 
findings reported in Reference [7]. These results highlight 
the importance of considering the zero-point energy and the 
neck-breaking probability parameter in the three-dimensional 
collective model for a more accurate description of the frag- 
ment mass distribution in fission processes. 


VII. CONCLUSION 


In summary, this article presents a comprehensive analysis 
of the fission process in Th, U, Pu, and Cm isotopes using a 
Yukawa-Folded mean-field plus standard pairing model. By 
employing a macroscopic-microscopic framework, the poten- 
tial energy surfaces, fission paths, barriers, and fragment mass 
distributions were calculated. The study focused on investi- 
gating the impact of the pairing interaction on the mass dis- 
tribution of fission fragments. 
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Fig. 10. (Color online) Mass yields of 7“°Pu with different values of 
the zero-point energy Eo in Eq. (15), and neck-breaking probability 
parameter d in Eq. (18). 


Our results demonstrate that the pairing interaction plays 
a crucial role in shaping the fission process or ?30Th, 734U, 
240Pu, and 246Cm. The strength of the pairing interaction was 
determined by fitting experimental data of odd-even mass dif- 
ferences and barrier heights, which led to a better agreement 
between theory and experiment. Furthermore, we found that 
the fission fragment mass distribution is highly sensitive to 
changes in the pairing interaction strength for ?30Th, 734U, 
240Pu, and ?46Cm, with stronger pairing interactions favor- 
ing symmetric fission and weaker interactions leading to more 
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asymmetric fission. The odd-even mass differences for ?°°Th, 
2340, 240Pu, and ?46Cm at symmetric and asymmetric fission 
points were compared with experimental values, providing 
additional support for the findings regarding the role of the 
pairing interaction. 


Moreover, by comparing our theoretical calculations with 
experimental data, we have confirmed the accuracy of our 
model in describing the fission fragment mass distributions 
for Th, U, Pu, and Cm isotopes. The peak heights, widths, 
and positions of the fragment mass distributions are well- 
reproduced, demonstrating the effectiveness of our approach. 


Additionally, the study explores the effects of the zero- 
point energy parameter and neck-breaking probability pa- 
rameter on the fragment mass distribution for ?40Pu. It is 
observed that the zero-point energy primarily influences the 
peak value of fission fragments, while the neck-breaking 
probability parameter affects the position of the asymmetric 
peak. 


In conclusion, this research contributes to our understand- 
ing of the fission process by emphasizing the crucial role of 
the pairing interaction and its relationship with nuclear elon- 
gation. The well agreement between theoretical calculations 
and experimental data, along with the analysis of additional 
parameters, strengthens the validity and applicability of the 
proposed model. The insights gained from this study can in- 
form future investigations in the field of nuclear fission, ulti- 
mately advancing our knowledge of this fundamental process. 
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